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A wave function picks up, in addition to the dynamic phase, the geometric (Berry) phase when 
traversing adiabatically a closed cycle in parameter space. We develop a general multidimensional 
theory of the geometric phase for (double) cycles around exceptional degeneracies in non-Hermitian 
Hamiltonians. We show that the geometric phase is exactly tt for symmetric complex Hamiltonians 
of arbitrary dimension and for nonsymmetric non-Hermitian Hamiltonians of dimension 2. For 
nonsymmetric non-Hermitian Hamiltonians of higher dimension, the geometric phase tends to tt 
for small cycles and changes as the cycle size and shape are varied. We find explicitly the leading 
asymptotic term of this dependence, and describe it in terms of interaction of different energy levels. 

PACS numbers: 03.65.Vf 



Non-Hermitian dissipative terms enter a quantum sys- 
tem Hamiltonian when studying non-isolated systems, 
e.g., efFective Hamiltonians describing decay of unsta- 
ble states. It turned out that the non-Hermitian physics 
differs dramatically from the Hermitian physics in the 
presence of degeneracies (energy level crossings), even 
if the non-Hermitian system is close to the Hermitian 
one 0, The most important degeneracy intrinsic 
to non-Hermitian Hamiltonians is the exceptional point 
(EP), at which two eigenvalues and corresponding eigen- 
vectors coalesce, as opposed to the diabolic point (DP) 
degeneracy of Hermitian operators, at which the eigenval- 
ues coalesce while the eigenvectors remain different. EP 
degeneracies have been observed in laser induced ioniza- 
tion of atoms 0, microwave cavities I3,l5l|, in "crystals of 
light" S]. Similar phenomena (where the Hamiltonian is 
substituted by a different system operator) are encoun- 
tered in optics of absorptive media [3 , acoustics B , elec- 
tronic circuits 0, and mechanical systems [Tolllll |. 

A wave function of a quantum system, whose param- 
eters undergo adiabatic cyclic evolution, acquires a com- 
plex factor dependent only on the loop in parameter 
space and, thus, called geometric or Berry phase [l^. 
Geometric phases in non-Hermitian systems were stud- 
ied in 0, H [ll [ll El El [Hill 13 . In such systems, it 
is important whether EP is inside the closed path or not. 
For Hamiltonians given by specific 2x2 matrices, the ge- 
ometric phase for a (double) cycle around EP was found 
to be exactly tt. Later this result was verified experimen- 
tally in 0.' So far, EPs have been observed in decaying 
systems described by symmetric effective Hamiltonians. 
This is the case when the corresponding isolated system 
is time-reversal (described by a real symmetric Hamil- 
tonian). Time-irreversal interactions, e.g., with exter- 
nal magnetic field, break the symmetry of the efFective 
Hamiltonian. 
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We should note that the existing theoretical studies 
for the geometric phase around EPs rely on the possibil- 
ity of reducing the system to the two-dimensional form. 
However, one should be aware that the geometric phase 
is generally not preserved under such reduction, as this 
reduction is given by a parameter dependent change of 
basis. For example, we mention the change of geomet- 
ric phase under the parameter-dependent magnetic gauge 
transformation [l3 |. 

In this paper, we develop a general multidimensional 
theory for geometric phases around EPs. We show that, 
for symmetric complex Hamiltonians of arbitrary dimen- 
sion and for general non-Hermitian Hamiltonians of di- 
mension 2, the geometric phase is exactly tt. However, 
for nonsymmetric non-Hermitian Hamiltonians of higher 
dimension, the geometric phase generally diverges from 
TT as the cycle size increases. We find explicitly the lead- 
ing term of this divergence. It describes the background 
influence of energy levels not involved in the EP degener- 
acy. We note that the divergence from tt is related to irre- 
versible Hermitian terms, rather than to non-Hermitian 
dissipative terms. 

Let H{X) be a non-Hermitian complex Hamiltonian 
smoothly dependent on a vector of m real parameters 
X = {Xi, . . . , Xm)- For simplicity, we consider Hamilto- 
nians represented by non-Hermitian complex matrices of 
arbitrary dimension, but the results are valid in infinite 
dimensional case as well. Let En{X) be the eigenvalues 
of H{X) (labeled n), and \'ipn{X)) be the corresponding 
eigenvectors. In multiparameter space, a set of EPs de- 
fines a smooth surface of codimension 2 j23 • For clarity, 
we assume that the number of parameters is three (then 
EPs form a curve), keeping in mind that the results be- 
low are valid for any number of parameters. Consider the 
EP curve, corresponding to the coincidence of the levels 
En = En+1 and the eigenvectors \ipn{X)) = |V'n-i-i(-'^))- 
Let C = {X{t) : 0<i<T}bea cycle making one turn 
around this EP curve in parameter space, see Fig.^ We 
assume that there are no degeneracies (multiple eigen- 
values) at points of the cycle C, as well as there are 
no other EP curves inside C. We note that EP is the 
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FIG. 1: Cycle around EP in parameter space. 



only generic codimension 2 degeneracy for complex non- 
Hermitian Hamiltonians smoothly dependent on param- 
eters j2(]|. Thus, strictly speaking, EP is the only degen- 
eracy that can be encircled by C in generic systems. 

Let |«'„(0)) = \MXm) and |«'„+i(0)) = 
\%pn+i{X (Q))) be the interacting quantum states aX t — 
0. After traversing the cycle C once, the states in- 
terchange (up to the phase multiplier) 0. When 
making two turns, both states return to their initial 
values picking up, in addition to the usual dynamical 
phase 5n = —ji J^^ E„{t)dt, a geometric phase 7„ [T^ : 
|*„,„+i(2T)) = e*(*'"+^")|*„,„+i(0)). Note that, due to 
the interchanging of the states, we have Sn — (5„+i and 
7n = 7ra+i- For non-Hermitian systems, the geometric 
phase is given by the integral [Tj, [Tjl 



7n = 7n+l 



{Mx)\dMx)) 
2c {Mx)\Mx)) ' 



(1) 



where {ipn{X)\ is the left eigenvector corresponding to 
En{X). The integral in ^ is evaluated over the cycle 
C traversed twice in the increasing time direction (we 
denote this by 2(7). The right and left eigenvectors are 
orthogonal at EP [21| , which means that the denominator 
of the integral expression in ^ is zero at EP. 

First, consider complex symmetric Hamiltonians: 
H{X) = H^{X). In this case the left and right eigenvec- 
tors are complex conjugate: (V'„,ti+i(-'^)| = {'4'n,n+i{X)\. 
By using this property, we write in the form 



In = I 



{Mx)\dMx)) 

2C (M^MniX)) 



(2) 



d\n{^P„{X)\MX)). 
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The phase 7„ is equal to the change of the complex quan- 
tity In {il;n{X)\ijjn{X)) over the double cycle 2C. This 
change depends only on t he number of turns made by the 
complex number (?/;„(X)|-0„(X)) around zero in complex 
plane, where zero is a branch point of the complex loga- 
rithm function. In one turn around zero, the logarithm 
changes by ±27ri p^ . 

Since the geometric phase 7„ does not depend on the 
form and size of the cycle, we evaluate 7„ by considering 



small cycles C around a point X-^-p of the EP curve. At 
Xep, two eigenvalues of the Hamiltonian iJ^p = H{Xj^p) 
coalesce: Eep — £'„(Xep) — £'„+i(A'ep). Eep has a sin- 
gle eigenvector |xo'') = |?/'„(Xep)) = |V'„,4-i (X^p)) and an 
associated vector \xi^) determined by [ifll 

HEp\xr)=EEp\xT), HEp\xT)^EEp\xr) + \xT)- (3) 



The left eigenvector (xo''| = (V'n(-^Ep)l = (V^„+i(Xep)| 
and associated vector (xf | are determined by 



{xT\HEP^EEp{xri {xT\HEP^EEp{xT\ + {xri (4) 

Recall that (Xo'^IXo^) = 0- Additionally, we impose the 
normalization conditions 

(xrm = (xgnxr) = i, (xrixr) = o. (5) 

In the neighborhood of X^p, we have 



\^n,^+lix)) - \xr)±vJ^\xr)+oW\\x - x^ 



(V„,„+i(X)| = (xg^|±VM(xri+o(VII^-^Bp||), 
where fi is the linear scalar function of parameters 



(6) 



m(X) = J2 {xr\dH/dX,\xr){Xj-X 
i=i 



EP\ 



(7) 



with the derivatives taken at ATep; the equation fi = 
gives the tangent of the EP curve in parameter space . 
By using jSJ), 10, and the property (Xo^il = (Xo^il for 
symmetric matrices, we obtain 



{MX)\MX)) = 2^+o{^\\X -XepW). (8) 

The complex number fj, makes one turn around zero in 
complex plane for one cycle C in parameter space. Hence, 
{ipn{X)\ipn{X)) makes a single closed loop around zero in 
complex plane for the double cycle 2C. As a result, the 
complex logarithm function in jSJ changes by ±27ri, and 
we obtain 7„ — ±tt. The sign depends on the direction of 
the cycle in complex plane; it does not influence the final 
result, since the phase is determined up to the additional 
term 27rfc for any integer k. 

There is the geometric phase analogy between EPs of 
complex symmetric Hamiltonians and DPs of real sym- 
metric Hamiltonians. For real symmetric Hamiltonians, 
just like for complex symmetric Hamiltonians, the geo- 
metric phase is "produced" only by the degeneracies: it is 
TT if the degeneracy is encircled, and zero otherwise |14| . 
Such phases, which do not depend on the shape (geome- 
try) of the cycle, are called topological [23|. The major 
difference between complex and real cases is that the cy- 
cle should be traversed twice for EP and once for DP. 
When a complex symmetric perturbation is given to a 
real symmetric Hamiltonian, DP splits into two EPs |2^ . 
One can say that each EP takes half of the geometric 
phase of DP (counted per single cycle). 
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Now, consider nonsymmetric non-Hermitian Hamilto- 
nians. We study the local structure of the EP degener- 
acy by means of the versal deformation theory of matri- 
ces [13,12^. The eigenvectors \i'n[X)) and \ipn+i{X)) 
are nonsmooth functions of parameters at X-^-p. How- 
ever, together they define a two-dimensional invariant 
linear subspace, which smoothly depends on parameters. 
This invariant linear subspace can be given by two vec- 
tors |xo(-^)) and smoothly dependent on param- 
eters: |xo,i(-^)) are linear combinations of \4'n{X)) and 
\il}n+i{X)) and satisfy the equations 26] 



H{X)\xa{X)) = s{X)\xa{X)) +p{X)\xi{X)), 
H{X)\xi{X)) - s{X)\xi{X)) + \xo{X)). 



(9) 



Here s{X) = {En{X) + En+i{X))/2 and p{X) = 
{En+i{X) ~ En{X))'^ / 4 are smooth scalar functions. At 



X = X-Ep, where s{Xep) = E^p and p{X-e 



0, ® 



yield the Jordan chain equations |j2Jl. Hence, Ixo(-'^ep)) 
IXo^) is the eigenvector and \xi{Xj^p)) = |xi^) is the as- 
sociated vector of the double eigenvalue E^p ■ By means 
of ©, the eigenvalues En^n+i{X) and corresponding 
eigenvectors are found as 



En^n+l{X) 
\^n,n+l{X)) 



six) ± 



\XoiX))±y/piX)\xiiX)), 



(10) 



where two Riemann sheets of the complex square root 
correspond to En{X) and En+i{X). We remark that the 
function ii{X) in ||7J) is the linearization p{X) at X^p- 
Similarly, the vectors (xo,i(^)l are introduced for the left 
eigenspace: they determine the left eigenvectors as 



(^n,n+i(^)l = (Xo(^)l ± yMX){xo{X)\, (11) 
and satisfy the orthonormality conditions 

(Xo(^)IXo(^)) = (xi(^)IXi(^)) =0, 



(Xi(X)|xo(^)) = (Xo(^)IXi(^)>=l. 



(12) 



At EP, (xo^l = (xo(-'^Ep)l is the left eigenvector and 
(Xi^l — (xi(A'ep)I is the left associated vector. 
By using 1^01)- (Q^l in CJ, we obtain 



In = 77 



2C 



+ 1 



2C 



{XoiX)\dxo{X))+p{X)(xi{X)\dxi{X)) 
2VpiX) 



(13) 



+ 2 f {{xoiX)\dxiiX)) + {xi{X)\dx,{X))). 

The double cycle 2C corresponds to a single cycle of the 
square root ^p{X) around zero in complex plane. Hence, 
the first integral in (|13|l equals ±27ri, where the sign de- 
pends on the direction of the cycle in complex plane. The 
second integral in (|13|l vanishes, since the square root in 
the denominator has opposite signs when traversing the 



first and second cycles. Finally, the third integral is the 
same for the first and second cycles. As a result, we have 



7„ = ivr + * y ((xo(X)Mxi(^))+(Xi(^)Mxo(^))). 

(14) 

Remark that the integral in (|14|l is taken over one cycle 
C in the increasing time direction. 

First, consider Hamiltonians given by 2 x 2 general 
complex matrices. According to (|12() . the 2x2 matrix 
|l)(xi(X)| + |2)(xo(X)| is the inverse of |xo(^))(l| + 
|xi(X))(2|, where |1) = (1, 0) and |2) = (0, 1) are 
the unit vectors. Hence, components of the vec- 
tors (xo,i(Ar)| can be expressed explicitly in terms of 
the components of |xo,i(A^))- By using these expres- 
sions, we transform the integral in H14|l to the form 
dlndet(|xo(-'^))(l| + lxi(-'^)) (2|); it vanishes since the 
2x2 matrix |xo(Ar))(l| -I- |xi(-'^))(2| is everywhere non- 
singular by definition. Hence, for 2x2 general non- 
Hermitian Hamiltonians, the geometric phase equals ±7r 
and does not depend on the loop shape, similar to the 
case of symmetric complex Hamiltonians. This result 
justifies the existence of topological indices describing the 
polarization ellipses around C points in crystal optics 0] . 

For multidimensional non-Hermitian Hamiltonians, 
the integral in (I14II is generally nonzero. Consider a cy- 
cle C = {X{t) = Xep + eX{t) : < i < T } making one 
turn around EP, where £ is a small positive parameter 
controlling size of the cycle. Formulae for derivatives of 
|Xo,i(^)) and (xo,i(A')| at X^p are provided by the ver- 
sal deformation method |26l | . By using these formulae in 
()14|l . we obtain the asymptotic expression 



In 



±7r- 



lae 



0(e3 



(15) 



where the complex constant a is given by the integral 

a= <f (2{xr\HiiG-'~\xT){xT\)dH,\xr) 

Jc ^ (16) 

+ {xr\HiG-^dHi\xr)+{xr\HiG~'dH,\xr) 



Here Hi{X) = Y.%i{dH / dXj)Xj and dHi{X) = 
'Yl^^iidH / dXj)dXj with the partial derivatives taken at 

Xep, and G — H^p — Ej^pl + |xi^)(Xi''l is a nonsingular 
matrix (/ is the identity operator). The correction term 
iae^ is determined by the information about the system 
at EP (this includes eigenvectors, associated vectors, and 
first derivatives of the Hamiltonian with respect to pa- 
rameters) and by the cycle shape X{t). Details of the 
derivation of (|15|l will appear elsewhere |23 . 

The physical meaning of the constant IjlGfl can be un- 
derstood by using the eigenvector expansion of the unity 
and of the Hamiltonian at EP: 



/ = ixn(xf i + ixrxxr 
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= \xr){xr\+E^A\xr){xr\+\xr){xr\) 



E 



(18) 



where i;^^ = Sfe(XEp), IVD = l^fc(^Ep)), and (^-^1 = 
('(/'fe (Xep ) I . Here we assume the normahzation condition 
for the left and right eigenvectors (V'fc^'IV'fe^) = 1- Recall 

that m^xZ) = (xg'il^D = and (V^ff IV'D = if 
k =/= k'. Expression 1)18(1 represents the transformation 
of Hep to the canonical Jordan form written in terms of 
eigenvectors and associated vectors ,21] . By substituting 
(|17|l and (|18|l into the expression for the matrix G, after 
a series of manipulations, we transform pt)|) to 



{xr\Himn{i^T\dHi\xf 




{El- - EepY 
{xT\HMT)m''\dHi\xf 



Ee 



(19) 



{xf\HMT)m''\dH,\xf 



(E, 



EP 

k 



E^, 



The terms (xoril^il^D and {^r\dHi\xfi^i) describe 
the interaction of the degenerate level E^p with the lev- 
els ii^fe, k ^ n,n + 1 at the EP. Thus, the change of the 
geometric phase with the cycle size and shape variation 
is due to the influence of the energy levels not involved 
in the EP degeneracy. One can see that, if the difference 
Er - E EP is big, the influence of the level Ek is pro- 



portional to {Ef - — _Eep)~^ and can be neglected. How- 
ever, if Ef- — Eep is small, the change of the geometric 
phase due to the interaction with Ek grows proportion- 
ally to {Ef- — Eep)~^ and may be big. In the extreme 
case Ef- — E^p 0, i.e., near the triple degeneracy 
En = En+i = Ek, we have a oo. Hence, triple degen- 
eracies require special investigation. 

Thus, for nonsymmetric non-Hermitian Hamiltonians, 
the deviation of the geometric phase from tt is the mul- 
tidimensional phenomenon, which cannot be captured in 
two-dimensional approximations. Asymptotic expression 
p5|l with the coefficient ((16|l for the correction term was 
confirmed by numerical simulations for particular Hamil- 
tonians of dimensions 3 and 4. We believe that this 
change of the geometric phase, which is intrinsic to non- 
symmetric non-Hermitian Hamiltonians, can be verified 
in future experiments. For example, this effect should ex- 
ist when studying the decay of nearly degenerate unsta- 
ble states for time-irreversal systems (Hamiltonians must 
have nonsymmetric Hermitian terms). Probably, the ex- 
perimental approach of can be used for this purpose 
if one manages to break the symmetry of the effective 
Hamiltonian in a controllable way. 

Based on the expansions of eigenvectors near EP we 
have shown that in the general case the geometric phase 
integral can be evaluated by methods of complex analy- 
sis. This can be regarded as a response to Arnold |2a| 
who suggested to develop a theory of "residues" to cal- 
culate the Berry phase. 
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